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In  this  study,  a  full-infinite  interval-stochastic  mixed-integer  programming  (FIMP)  method  is  developed 
for  planning  carbon  emission  trading  (CET)  under  dual  uncertainties.  FIMP  has  advantages  in 
uncertainty  reflection  and  policy  analysis,  particularly  when  the  input  parameters  are  provided  as 
crisp  and  functional  intervals  as  well  as  probabilistic  distributions.  The  developed  FIMP  is  applied  to  a 
real  case  study  for  managing  carbon  dioxide  (C02)  emissions  with  trading  scheme  of  Beijing's  electric 
power  system  (EPS).  Electric  power  industry  is  one  of  the  major  sources  of  C02  emission  in  China.  It  is 
essential  to  accumulate  relevant  experience  to  provide  a  reliable  basis  for  establishing  a  regional  or 
national  CET  market,  so  as  to  prepare  for  docking  with  the  international  market.  This  is  the  first  attempt 
to  introduce  CET  scheme  into  Beijing's  EPS  to  mitigate  C02  emissions.  The  solutions  for  energy  supply, 
electricity  generation,  carbon-quota  allocation,  and  capacity  expansion  are  obtained.  They  cannot  only 
be  used  for  formulating  C02-reduction  policies  and  assessing  the  associated  economic  implications  in 
purchasing  emission  permits  or  bearing  economic  penalties,  but  also  facilitate  analyzing  various 
policies  when  pre-regulated  electricity-generation  plans  and  pre-defined  C02-emission  schemes  are 
violated. 
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1.  Introduction 

Damage  caused  by  global  warming  is  happening  far  faster  than 
experts  predicted  or  anticipated.  Potential  threats  generated  by 
global  warming  may  contain  the  increase  in  surface  temperature, 
the  change  in  the  global  climate,  the  rise  in  ocean  level,  even  the 
disruption  in  food  production  [1,2].  The  preponderance  of  emer¬ 
ging  scientific  evidence  emissions  of  heat-trapping  carbon  dioxide 
(C02)  from  fossil  fuels  are  tracking  near  the  highest  level  con¬ 
sidered  by  the  Intergovernmental  Panel  on  climate  change  [3], 
Generally,  C02  emission  is  closely  related  to  human  activities 
such  as  power  generation,  deforestation,  transportation,  as  well 
as  industrial,  residential  and  commercial  activities.  C02  is  dis¬ 
charged  into  the  atmosphere  from  combustions  of  fossil  fuels 
such  as  oil,  natural  gas,  and  coal  as  sources  of  energy.  Since  the 
Industrial  Revolution,  human  activity  has  increased  the  amount  of 
greenhouse  gas  (GHG)  in  the  atmosphere,  leading  to  increased 
radiative  forcing  from  GHG.  Amongst  several  environmental 
pollutants  causing  climate  change,  C02  is  responsible  for  58.8% 
of  the  GHG  [4].  Meanwhile,  the  concentration  of  C02  has 
increased  by  36%  since  1750  [5].  Specially,  fossil  fuel  burning  in 
electric  power  system  (EPS)  has  produced  about  three-quarters  of 
the  increase  in  C02  from  human  activity  over  the  past  20  years, 
which  is  one  of  the  major  contributory  factors  for  occurring  global 
warming  [3],  Since  it  is  undeniable  that  many  countries  are  still 
mainly  dependent  on  fossil  fuels  to  sustain  their  energy  genera¬ 
tion  and  supply,  C02  emission  will  exacerbate  global  warming 
and  lead  to  environmental  destruction  and  health  hazards  which 
are  already  quite  rampant  nowadays. 

Due  to  the  harmful  impact  of  excessive  C02  emission,  regula¬ 
tory  agencies  have  enacted  strict  regulations  to  limit  their  emis¬ 
sions.  As  a  crucial  international  agreement,  Kyoto  Protocol  was 
signed  in  1997  to  respond  to  global  warming  through  reducing 
GHG  emission.  Recently,  a  milestone  is  marked  in  the  Kyoto 
Protocol  with  the  2009  Climate  Summit  in  Copenhagen,  Denmark, 
with  all  participating  countries  further  committed  themselves 
in  fulfilling  the  protocol’s  obligations  before  the  commitment 
period  due  in  2012  [6].  Currently,  there  are  several  CET  markets 
throughout  the  world,  including  clean  development  mechanism 
(CDM)  project  market,  European  Union  greenhouse  gas  emission 
trading  scheme  (EU  ETS),  UK  emissions  trading  group  (ETC), 
Chicago  climate  exchange  (CCX),  and  Australian  climate  exchange 
(ACX).  In  2005,  EU  ETS  was  proposed  as  an  allowance-based 
transaction  which  referred  to  the  excess  emission  reduction 
trading  under  the  total  amount  control  among  countries  identi¬ 
fied  by  Kyoto  Protocol.  As  a  major  pillar  of  EU  climate  policy,  EU 
ETS  is  the  first  large  emission  trading  scheme  in  the  world  [7].  In 
2010,  emission  trading  revenue  of  EU  ETS  reached  to  $  119.8 
billion,  accounting  for  84%  of  the  income  of  the  global  carbon 
trading.  EU  ETS  covers  11,000  power  stations  and  industrial 
plants  in  30  countries,  which  are  collectively  responsible  for 
40%  of  its  total  GHG  emissions.  Under  EU  ETS,  all  members  agree 
on  their  national  emission  caps  which  have  to  be  approved  by  the 
EU  commission,  and  then  allocate  allowances  to  their  industrial 
operators,  and  validate  the  actual  emissions  in  accordance  with 
the  relevant  assigned  amount  [8]. 

As  a  market-based  scheme,  emissions  trading  could  lead  to  a 
desired  market  transformation  and  encourage  the  application  of 


energy  efficient  technologies,  as  well  as  ensure  the  achievement 
of  emission  reduction  target  [9].  During  the  period  of  2005-2009, 
the  EU  ETS  successfully  reduced  carbon  emissions  up  to  5%,  with  a 
limited  economic  impact  of  less  than  1%  of  total  gross  domestic 
product  (GDP).  The  major  British  industries  on  average  faced  a 
cost  increase  amounting  to  4%  at  a  carbon  price  of  $28.8  per 
metric  ton  [10].  Although  no  federal  emissions  trading  market 
in  United  States,  regional  markets  are  in  place  or  are  under 
development,  including  the  Regional  Greenhouse  Gas  Initiative, 
Western  Climate  Initiative,  and  the  Midwestern  Greenhouse 
Gas  Accord.  According  to  these  fundamental  policies,  a  number 
of  approaches  are  available  to  curb  emissions  and  tackle  climate 
change.  For  example,  the  cap  and  dividend  system  proposed  by 
Peter  Barnes  is  a  market-based  approach  to  reduce  carbon 
emissions  without  reducing  household  incomes.  Under  this 
mechanism,  fossil  fuel  supplies  and  emissions  are  capped.  Pollu¬ 
ters  have  to  pay  a  premium  to  emit,  which  can  be  returned  back 
to  all  citizens  equally  [11].  In  another  alternative  approach, 
polluters  voluntarily  exchange  the  emissions  within  the  cap  and 
trade  mechanism  and  pay  extra  taxes  when  emissions  overrun 
the  cap.  Existing  research  works  mainly  focused  on  economic 
input-output  modeling,  life  cycle  assessment,  and  operating 
research  models  to  examine  the  policy  implications  [12]. 

China  is  the  second  largest  GHG  emitter  throughout  the  world, 
which  pledges  to  reduce  emissions  per  unit  of  economic  output 
by  40-45%  relative  to  2005  level.  However,  with  rapid  economic 
development,  China’s  energy  consumption,  especially  electricity 
consumption,  has  grown  swiftly.  Coal-fired  power  occupies  more 
than  75%  of  the  total  electricity  generation  in  China,  which 
presents  particular  phenomenon  all  around  the  world.  Addition¬ 
ally,  the  amount  of  electricity  generated  from  coal-fired  power 
plant  increased  rapidly  over  recent  years.  For  example,  from  1991 
to  2005,  the  amount  has  been  maintained  at  90-96%,  which  means 
that  large  amounts  of  C02  have  been  emitted.  In  2005,  carbon 
emissions  from  the  Chinese  power  sector  have  reached  38.73%  of 
total  emissions  of  primary  energy  [13].  Electric  power  industry  is 
one  of  the  major  sources  of  C02  emission  in  China.  It  is  funda¬ 
mental  to  accumulate  relevant  experience  to  provide  a  reliable 
basis  for  establishing  regional  and/or  national  CET  markets  [14]. 
However,  Chinese  current  CET  market  is  a  CDM  one  [15],  Many 
cities  (e.g.,  Beijing,  Tianjin  and  Shanghai)  have  set  up  the  environ¬ 
ment  or  property  exchange  as  the  energy  exchange  trading 
platform  for  CDM  projects.  According  to  the  studies  of  World  Bank, 
China  has  potential  for  C02-emission  reduction  from  100  million 
ton  to  200  million  ton  per  year.  It  can  provide  more  than  half  of  the 
global  CDM  projects.  In  May  1998,  China  signed  the  Kyoto  Protocol 
at  United  Nations  Headquarters,  and  in  August  2002,  formally 
approved  the  Kyoto  Protocol,  which  means  that  China  start  opera¬ 
tion  of  the  CDM  projects.  To  enhance  the  effective  management  of 
the  CDM  projects  and  to  ensure  that  their  projects  go  into 
an  orderly  manner,  Chinese  government  promulgated  interim 
measures  for  the  operation  and  management  of  the  clean  devel¬ 
opment  mechanism  project  in  June  2004.  In  October  2005,  interim 
measures  was  amended,  clean  development  mechanism  project 
management  approach  is  enacted.  However,  current  CDM  market 
is  mainly  a  market  for  buyers.  As  a  market  for  sellers,  the  project 
owners  in  developing  countries  are  weak  of  negotiation  right, 
which  causes  them  being  disadvantaged.  At  the  same  time, 
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it  is  also  imperative  to  realize  energy  saving  and  C02-emission 
reduction,  respond  to  international  pressure  and  participate  in 
international  conventions  roundly.  Moreover,  few  research  works 
concerning  the  impacts  of  CET  on  EPS  of  China  have  been  proposed 
in  a  computable  framework  [14,16,17]. 

Previously,  various  mathematical  models  were  undertaken  for 
tackling  the  issues  of  CET;  some  of  them  were  introduced  for 
analyzing  issues  related  to  effectiveness  and  equality  in  CET 
[18-22].  For  example,  Bosello  and  Roson  [18]  formulated  an 
integrated  assessment  model  to  analyze  the  CET  and  equity  in 
international  agreements,  which  explored  the  distributional  con¬ 
sequences  of  alternative  emission  trading  schemes.  They  found 
that  the  introduction  of  a  competitive  market  for  emission 
permits,  especially  when  this  market  has  no  developing  countries, 
would  dramatically  lower  total  abatement  costs.  Holtsmark  and 
Maestad  [19]  employed  a  numerical  model  to  assess  the  signifi¬ 
cance  of  international  emission  trading  for  the  oil,  coal  and 
gas  markets.  Yi  et  al.  [20]  constructed  an  intensity  allocation 
model  based  on  three  indicators  (i.e.,  per  capita  GDP,  accumulated 
fossil  fuel  related  C02  emissions  and  energy  consumption  per  unit 
of  industrial  added  value),  where  equity  principles  for  target 
allocation  were  taken  into  account.  Lu  et  al.  [21]  presented  a 
carbon  regulation  based  duopoly  model  to  evaluate  the  effective¬ 
ness  and  equity  of  various  carbon  policies  including  emission 
standard,  carbon  tax,  and  emission  trading. 

Many  researchers  evaluated  the  economic  efficiency  of  redu¬ 
cing  carbon  emission  which  was  the  goals  of  the  EU  ETS  were 
criticized  for  not  being  ambitious  enough  [22],  Since  it  was 
generally  either  technically  infeasible  or  economically  impossible 
to  design  processes  leading  to  zero  emission  of  C02,  authorities 
and  decision  makers  always  sleeked  to  control  the  C02  emission 
to  level  at  which  the  effect  was  minimized  [23].  Therefore, 
it  is  imperative  to  make  a  criterion  of  allowable  levels  of  C02 
emissions.  Mathematical  models  were  improved  ceaselessly  to 
response  the  changing  factors  related  to  emission  permits 
[24-27].  For  example,  Haurie  and  Viguier  [24]  proposed  a  com¬ 
putable  stochastic  equilibrium  model  to  represent  possible 
competition  between  Russia  and  China  on  the  international 
market  of  carbon  emission  permits;  they  analyzed  the  impact  of 
this  competition  on  the  pricing  of  emission  permits  and  on  the 
effectiveness  of  the  Kyoto  and  post-Kyoto  agreements,  without 
the  United  States  participation.  Rehdanz  [25]  developed  a  two- 
country  game  model  to  analyze  the  coordination  of  domestic 
markets  for  tradable  emission  permits,  where  countries  deter¬ 
mined  their  own  emission  reduction  targets.  Bristow  et  al.  [26] 
developed  a  pooled  model  to  explore  the  influence  of  key 
design  attributes  on  the  acceptability  of  a  personal  carbon  trading 
scheme,  where  permit  prices  vary  among  different  models  because 
of  differences  in  percentage  reductions  in  emissions,  baseline 
prices  of  carbon-energy,  and  price  sensitivity  of  users  [27]. 
In  addition,  ICL1PS  (Integrated  Assessment  of  Climate  Protection 
Strategies)  Climate  Model  (ICM)  can  be  used  to  simulate  the 
atmospheric  retention  and  metabolism  of  C02.  Linking  the 
economic  and  the  climate  impact  components  in  the  ICLIPS 
framework,  major  greenhouse  gases  (including  C02,  CH4,  N20, 
halocarbons,  SF6,  tropospheric  and  stratospheric  03,  and  strato¬ 
spheric  water  vapor)  as  well  as  the  radiative  effects  of  aerosols 
originating  from  S02  emissions  and  from  biomass  burning  were 
considered.  Several  researchers  analyzed  the  equity  issues  and 
cost-effective  emission  trajectories  based  on  the  ICM  [28-30], 

Carbon  emissions  are  easier  to  monitor  and  control  from  a 
limited  number  of  centralised,  large  power  stations  than  from 
millions  of  vehicles,  small  boilers  and  even  ruminant  animals. 
Accordingly,  electricity  sector  has  become  a  prime  target  where 
carbon  emission  controls  are  implemented  and  carbon  emission 
mitigation  is  valued  [31].  Previously,  a  number  of  researchers 


managed  C02  emissions  of  EPS  through  CET  scheme  [32-37]. 
For  example,  Bonacina  and  Gulli'  [32]  analyzed  the  impact 
of  C02-emission  allowance  trading  on  electricity  pricing  in  short 
term.  Bernard  et  al.  [33]  proposed  a  computable  dynamic  game 
model  of  the  strategic  competition  between  Russia  and  develop¬ 
ing  countries  (mainly  represented  by  China).  Linares  et  al.  [34] 
presented  an  oligopolistic  generation-expansion  model  for  the 
electricity  sector.  Chappin  and  Dijkema  [35]  presented  an  agent- 
based  model  to  elucidate  the  effect  of  CET  on  the  decisions  of 
power  companies  in  an  oligopolistic  market;  in  power  generation, 
the  economic  effect  of  CET  was  not  sufficient  to  outweigh  the 
economic  incentives  in  choosing  coal.  Sadegheih  [36]  proposed  an 
optimization  model  to  search  for  solutions  to  power  network 
planning  under  the  CET  program,  which  possessed  the  ability  to 
minimize  the  total  cost  with  cost-effective  and  environmentally 
friendly  manner.  Koo  et  al.  [37]  proposed  a  robust  optimization 
approach  to  plan  sustainable  energy  systems  of  South  Korea, 
which  could  help  decision  makers  to  determine  the  optimal 
capacities  of  power  plants  and/or  carbon  capture  and  storage,  as 
well  as  volumes  of  emissions  trading  that  could  meet  the  required 
emission  level  and  satisfy  energy  demand  from  various  user- 
sections  with  minimum  costs  and  maximum  robustness. 

In  practical  CET  planning  problems,  a  variety  of  complexities 
and  uncertainties  exist  among  different  electricity-generation 
activities  and  their  socio-economic  and  environmental  impli¬ 
cations  [38-41].  In  CET  programs,  uncertainties  often  arise 
due  to  inconstant  commitments  or  changes  of  emission  trading 
regulatory  [42-47].  For  example,  C02  discharged  from  various 
electricity-generation  activities  (e.g.,  power  plant,  storage  hydro¬ 
electric  plant,  or  biomass  power  plant)  can  be  influenced  by  some 
uncertain  events  (e.g.,  emission  limitation),  which  may  fluctuate 
time  to  time.  Meanwhile,  errors  in  estimated  modeling  para¬ 
meters  (e.g.,  economic  penalty)  could  be  possible  sources  of 
uncertainties.  These  complexities  have  placed  many  EPS  manage¬ 
ment  problems  beyond  the  conventional  optimization  methods. 
Therefore,  it  is  indispensable  to  inject  more  and  more  momentum 
to  the  EPS  planning,  including  considerations  for  diversity  of 
energy  activities,  structure  of  electricity  generations,  variation 
of  system  conditions,  uncertainty  of  impact  factors,  dynamics 
of  capacity  expansion,  as  well  as  the  associated  environmental 
implication  [48].  As  a  result,  a  few  research  efforts  were  con¬ 
ducted  for  dealing  with  various  uncertainties  in  CET  programs 
[46,49,50].  For  example,  Monni  et  al.  [49]  estimated  uncertainties 
in  different  emissions  trading  schemes;  according  to  their  results, 
uncertainty  in  emissions  from  the  EU15  and  the  EU25  included  in 
the  first  phase  of  the  EU  emissions  trading  scheme  (2005-2007) 
was  +  3%  (at  95%  confidence  interval  relative  to  the  mean  value). 
Chen  et  al.  [50]  developed  a  two-stage  inexact-stochastic  pro¬ 
gramming  model  for  planning  GHG-emission  trading  within  the 
electrical-power  systems.  Li  et  al.  [46]  proposed  an  interval-fuzz  y 
two-stage  stochastic  programming  method  for  planning  CET 
under  uncertainty,  where  three  trading  schemes  were  considered 
based  on  different  trading  participants. 

The  most  previous  optimization  models  were  formulated 
through  conventional  interval-parameter  programming  (IPP) 
approach.  Nevertheless,  the  conventional  IPP  has  difficulties 
when  the  right-hand  sides  of  a  model  are  highly  uncertain, 
especially  with  uncertainties  expressed  as  possibilistic  and/or 
probabilistic  distributions,  which  may  lead  to  the  loss  of  valuable 
information  in  many  real-world  decision-making  problems  [51]. 
Additionally,  IPP  can  only  solve  the  problems  containing  crisp 
interval  coefficients  [a,  b],  whose  lower-  and  upper-bounds  (i.e.,  a 
and  b )  are  both  deterministic  and  definitely  known.  Stochastic 
mathematical  programming  (SMP)  is  effective  for  decision  pro¬ 
blems  whose  coefficients  (input  data)  are  uncertain  but  could  be 
represented  as  chances  or  probabilities,  which  has  been 
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extensively  applied  to  energy  systems  planning  [52].  Two-stage 
stochastic  programming  (TSP)  is  a  typical  SMP  method,  which  is 
an  effective  alternative  for  tackling  problems  where  an  analysis  of 
policy  scenarios  is  desired  and  the  right-hand-side  coefficients  are 
random  with  known  probability  density  functions  (PDFs)  [53].  In 
practical  planning  problems,  it  is  important  to  expression  the 
relationship  between  mathematical  parameter  and  their  eco¬ 
nomic  influence  in  CET.  The  definition  of  crisp  interval  is  not 
suitable  for  all  cases  where  the  two  bounds  may  be  associated 
with  the  external  impact  factors  expressed  as  functional  interval 
[/(a),  g(i>)[,  where  /(a)  and  g(b)  are  the  functions  of  a  and  b, 
respectively.  An  attractive  technique  that  could  tackle  functional 
intervals  was  full-infinite  programming  (FIP),  which  was  pro¬ 
posed  to  deal  with  the  uncertainties  expressed  as  both  crisp 
intervals  and  functional  intervals. 

Therefore,  the  objective  of  this  study  is  to  develop  a  full-infinite 
interval-stochastic  mixed-integer  programming  (FIMP)  method 
in  response  to  the  above  challenges.  Then,  a  FIMP-based  carbon 
emission  trading  (FIMP-CET)  model  will  be  formulated  for  managing 
C02  emissions  of  Beijing’s  electric  power  system  (EPS)  with  trading 
scheme.  In  the  FIMP-CET  model,  emission-permit  limits  will  be 
allocated  to  each  power  plant  and  the  total  C02-emission  amount 
cannot  be  exceeded.  If  actual  emission  amount  less  than  the  allowed 
permit,  the  power  plant  can  sold  its  excess  C02  permit  to  achieve 
economic  benefit;  on  the  contrary,  it  has  to  buy  the  emission 
permits  from  the  others  or  subject  to  heavy  economic  penalties 
from  the  local  government.  The  modeling  results  will  be  used  for 
generating  a  range  of  decision  alternatives  under  various  system 
conditions,  and  thus  helping  decision  makers  to  discern  optimal 
power-generation  patterns,  and  gain  deep  insights  into  the  tradeoffs 
between  C02  emission  and  economic  objective. 


2.  Methodology 


Min  f  =  CT,X+J2phDT2Y  (2a) 

/i  =  i 

subject  to: 

ArX<Br,  r=l,2,...,m1  (2b) 

A[X+A'tY  >  wh,  teM,  M=  l,2,...,zn2,  h  =  \,2,...,v  (2c) 

Xj  >  0,  XjeX,  j=  1,2,. ...Hi  (2d) 

y,h  >  0,  yjheY,  j=  1,2,..., n2  (2e) 


Obviously,  model  (2)  can  deal  with  uncertainties  in  the  right- 
hand  sides  presented  as  probability  distributions  when  coeffi¬ 
cients  in  the  left-hand  sides  and  in  the  objective  function  are 
deterministic.  However,  in  real-world  decision-making  problems, 
the  quality  of  information  that  can  be  obtained  is  mostly  not 
satisfactory  enough  to  be  presented  as  probabilities.  Such  com¬ 
plexities  cannot  be  solved  through  model  (2).  Generally,  interval- 
parameter  programming  (IPP)  approach  is  effective  in  tackling 
uncertainties  expressed  as  interval  values  with  known  lower  and 
upper  bounds  but  unknown  distribution  functions.  Therefore, 
through  incorporating  IPP  and  TSP  within  a  general  optimization 
framework,  a  hybrid  two-stage  stochastic  programming  linear 
model  can  be  formulated  as  follows  [56]: 

Min/±  =C±X±+  J2PhDf2Y±  (3a) 

h  =  1 

subject  to: 

Ar±X±  kB*,  r=  1,  2,  ...,  m,  (3b) 


Two-stage  stochastic  programming  (TSP)  method  is  effective 
for  problems  where  an  analysis  of  policy  scenarios  is  desired  and 
the  related  data  are  mostly  uncertain  [54-56].  In  TSP,  the  first- 
stage  decision  is  to  be  made  before  uncertain  information  is 
revealed,  whereas  the  second-stage  one  (recourse)  is  to  adapt  to 
the  previous  decision  based  on  the  further  information;  the 
second-stage  decision  is  used  to  minimize  “penalties”  that  may 
appear  due  to  any  infeasibility  [47,50].  The  main  advantage  of  TSP 
is  its  capacity  in  dealing  with  recourse,  where  corrective  actions 
can  be  taken  after  a  random  event  has  taken  place.  A  general  TSP 
model  can  be  formulated  as  follows: 

z=  min  CTX+Em e a[Q(X,a))] 

subject  to: 

xeX  (la) 

with 

Q(x,co)  =  min /(m)7)/ 
subject  to: 

D(co)y>h(co)+T(co)x  (lb) 


where  Xek"1,  CeR "2  and  YeR "2.  Here,  co  is  a  random  variable 
from  space  ( Q,F,P )  with  Q  e  Rk,f  :  £2->R"2  ,h  :  £?->Rm2  ,D  :  C2-r 
Rm2  xn2 ,  an(j  T\Q-*Rm2xn  '.  By  letting  random  variables 
(i.e.,  co)  take  discrete  values  coh  with  probability  levels  ph 
(h  =  1,2,...,  v  and  Jjp h  =  1),  the  above  TSP  can  be  equivalently 
formulated  as  a  linear  programming  model  as  follows: 


At±X±  +At±Y±  <w±,  t=  1,2, ...,m2;  h  =  1,2 iz 


(3c) 


>0, 


eXd 


j=  1,2 . nt 


(3d) 


Yjt 


yJ^eY±,  j=  l,2,...,n2;  h=  1,2 . v 


(3e) 


where  A±  e  [R±  )m'xn\  x±  e(R±}m2X"2 


[R] 


1  xri) 


Dr  e  {R} 


1  xn2 


X±  e  (R  * } 


±  ,/!,  xl 


b; 


ElS1)"””1,  Cr±  6 
y±6{R±)"2Xl  and  {R*} 


denote  a  set  of  crisp  interval  parameters  and/or  variables;  super¬ 
scripts  ‘  -  ’  and  ‘  +  ’  represent  lower  and  upper  bounds  of  the  interval 
values,  respectively  [57-60].  In  model  (3),  the  decision  variables  are 
divided  into  two  subsets:  the  first-stage  decision  variables  X±  that 
must  be  determined  before  the  random  variables  are  disclosed,  and 
recourse  variables  Y  1  that  can  be  determined  after  the  random 
variables  are  disclosed.  In  the  model,  decision  variables  are  divided 
into  two  subsets:  Those  that  must  be  determined  before  random 
variables  are  disclosed  and  those  (recourse  variables)  that  will  be 
determined  after  the  uncertainties  are  disclosed. 

Model  (3)  can  handle  uncertainties  expressed  as  probability 
distributions  as  well  as  account  for  economic  penalties  with 
recourse  against  any  infeasibility;  however,  it  cannot  address 
uncertainty  expressed  as  functional  intervals  in  the  objective  and 
constraints.  In  real-world  decision-making  problems,  functional 
interval  is  defined  as  an  interval  with  its  lower  and  upper  bounds 
being  functions  of  an  independent  variable.  It  is  indicated  that  the 
lower-  and  upper-bounds  of  the  crisp  interval  can  vary  with  its 
independent  variables.  In  this  study,  crisp  interval  coefficients 
mean  the  conventional  intervals,  such  as  [a,  b],  whose  lower  and 
upper  bounds  (i.e.,  a  and  b)  are  both  deterministic  and  definitely 
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known.  This  is  based  on  the  assumption  that  these  interval 
coefficients  are  unchanged  even  if  they  could  be  affected  by 
associated  impact  factors.  Thus,  despite  the  effectiveness  of  the 
previous  methods  in  solving  models  containing  interval-para¬ 
meters,  they  are  still  limited  for  tackling  functional  intervals.  An 
attractive  technique  that  can  tackle  both  crisp  intervals  and 
functional  intervals  is  semi-infinite  programming  (SIP).  However, 
SIP  can  only  resolve  problems  with  merely  infinite  constraints  in 
deterministic  environment,  while  FIP  can  address  crisp  and 
functional  intervals  with  infinite  objectives  and  constraints.  Con¬ 
sequently,  when  the  coefficients  in  objective  function  and  con¬ 
straints  are  all  allowed  to  be  functional  intervals  and/or  crisp 
intervals,  a  full-infinite  interval-stochastic  mixed-integer  pro¬ 


gramming  (FIMP)  model  can  be  formulated  as  follows: 

V 

Min/1  =Cf/(T,)X±  +  E  PhDf12y±  for  all  t,-  e [zl,ru]  (4a) 

h  =  1 

subject  to: 

A1  (t,)X  1  <  B1  (Tf),  r=l,2,...,m1  (4b) 

At±X±+A't±Y±  <w±,  t=  1,2,. ..,m2;  h=  1,2,. ...v  (4c) 

x/*1  >  0,  xjteX±,  j=  l,2,...,ni  (4d) 

yjt  >  0,  yj^eY±,  j  =  1,2 . n2;  h=  1,2 . v  (4e) 


where  X  ±  (j=  1 ,  2 . n)  are  decision  variables;  C ^  (t,).  A1  (t,),  and 

B~  (Tj)are  functional  interval  parameters  in  objective  and  constraints. 
Cf  (Tj)  and  C/j(Ti)  (j=l,  2 are  positive  for  all  TiValues;  Cf  ( t,-) 

and  Cjf(T,)  (j=k+\,k+2 . n)  are  negative  functions  for  all  t,  values. 

As  a  result,  an  optimized  set  of  values  can  be  identified  by  having 
Ujt  as  decision  variables;  this  optimized  set  may  correspond  to 
minimized  system  cost  under  uncertain  generation  targets  x^  [56]. 
In  detail,  let  Xjf  =  xjt  +  AxjtUjt,  where  AxJt  =  Xj/  — xjj,  and  uJt  e  [0,1], 
When  xjf  approach  their  lower  bounds  (i.e.,  uJt  =  0),  a  relatively  low 
cost  would  be  obtained;  however,  a  higher  penalty  may  have  to  be 
paid  when  the  demand  is  not  satisfied.  Conversely,  when  Xj~  reach 
their  upper  bounds  (i.e.,  uJt  =  1 ),  a  higher  cost  would  be  generated 
but,  at  the  same  time,  a  lower  risk  of  violating  the  promised  targets 
(and  thus  lower  penalty).  Then,  model  (4)  can  be  transformed  into 
two  deterministic  submodels.  When  the  objective  is  to  be  minimized, 
the  submodel  corresponding  to  f~  can  be  firstly  formulated: 

k  i  n, 

Min/“  =  ^c/(Ti)x/  +  E  c/(Ti)x/ 
j=  1  j  =  ^1  +  1 

k2  v  k2  v 

+  EE Pmdi yjh+  E  E pjhdI yfh 

j  =l/]  =  l  j  =  k2  +  l/i  =  l 

for  all  t,  e  [ zl,ru ]  (5a) 

subject  to: 

k ,  n, 

E  I ad(Ti)  I +  Signfa^  (t,))x/  +  ^  | arJ( t,) |  ~Sign (a^tdlxT 

j  =  i  j  =  fci + 1 

<  br+  (t O.Vr  (5b) 


(f,  n, 

E  !  ay  |  +  Sign(  at+  )x/  +  ]T  |  arj  |  “  Sign(  )Xj+ 

j  =  i  j  =  fci  + 1 


+  E  I atj  I  Sign(a'+)yjh+  E  layl  sisn(a'tj  )y+ >  coh ,  Vt,h 

j  =  1  j  =  k2  +  l 


(5c) 

.X 
•  1 

IV 

o 

j=  1.2 . k, 

(5d) 

Xy+  >  0, 

j  =  ki  +  l,  k\  +2 . Hi 

(5e) 

yjh  ^  o, 

Vh;  j=l,2,...,k2 

(5f) 

y% 

V/t;  j  =  fe2  +  l,  k2  +  2,...,n2 

(5g) 

where  xA  J  =  1,2,...,  k,  are  interval  variables  with  positive  coeffi¬ 
cients  in  the  objective  function;  xjtj  =  k\  +  1,k\  +  2,...,ri\  are 
interval  variables  with  negative  coefficients;  yj~h>  0,  j=\,2,...,k2 
and  h=  1,2, ...,v  are  random  variables  with  positive  coefficients 
in  the  objective  function;  yjfc  >0,  j  =  fc2  +  l,fc2  + 2,...,n2  and  h  =  1, 
2,...,v  are  random  variables  with  negative  coefficients.  Solutions 

of  xjoPt(j=  I-2 . fei  )•  xUt(j=  h  +  1.  k ,  +2,...,n,),  y/tapt O' =  1.2 . k2) 

and  y+  O'  =  k2  +  l,fc2  +  2,...,n2)  can  be  obtained  through  submodel 
(5).  Based  on  the  above  solutions,  the  second  submodel  corre¬ 
sponding  to/+  can  be  formulated  as  follows: 

Min/+  =  Y^cY  (%i)x+  +  jr  c+( zt)x-  +  E  E  Piktfyjh  +  E  E  PM+>S 

j  =  1  j  =  fc,+l  j  =  l  h  =  1  j  =  k2+t  h  =  t 

for  all  t,  6  [t/,tu]  (6a) 

subject  to: 

k,  _  n, 

EKM  Sign(a^(T,))x/ +  J2  lQy(Ti)|  Sign(a t (t,))x/ <  b/(T,),  Vr 

j  =  i  j  =  fci  +  i 

(6b) 

ki  n , 

E  !  a‘J  I  s'Sn( ay  >XJ+  +  E  I  %  r  SiSn<  ay  >xjr 

j  =  i  j = fci + 1 


+  E  I ay'  I  siSn(fly  )y}£  +  E  I fl^r  Sign(a'+ )yj~h^co+ ,vt,h 

j  =  1  j  =  k2  + 1 


(6c) 

*j  ^xjoPt.  J  =  1,2,. ...ki 

(6d) 

0<x+  < x+pt,  j=  kr  + 1,  fci  +2, ...,ni 

(6e) 

yjh^yjhopr  v/i;  j = i , 2, .... /<2 

(6f) 

0<y+  <y+0pr  vh;  j  =  k2  + 1,  fc2+ 2 . n2 

(6g) 

solutions  of  x+pt(j  =  1,2, ...,fc, ),  xrpt0'  =  fci  +  l.fci  +2,...,^ ),  y+opt0’  = 
1,2,..., fci),  and  y/,opt,  Vb:  j  =  fc2  +  l,fc2  +  2,...,n2  can  be  obtained. 
Through  solving  submodels  (5)  and  (6),  interval  solutions 
associated  with  probability  information  can  be  obtained  as 
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follows: 

Xjopt  =  [Xjopt'Xjopti’^j 

(7a) 

y jhhopt  \yjhhoptJ 'jhhoptl 

(7b) 

/opt  =  \fopJopt]yj 

(7c) 

3.  Case  study 

3.2.  Statement  of  problems 

Beijing  is  the  capital  of  China;  with  the  increasing  economic 
development  and  population  growth,  the  amount  of  energy 
demand  has  rapid  growth  in  the  last  decades.  According  to 
Statistics  Bureau  of  Beijing,  total  energy  consumption  amounted 
to  42.29  million  tonne  of  coal  equivalent  in  2001,  while  65.70 
million  tonne  of  coal  equivalent  in  2009.  In  2010,  the  amount  of 
the  city’s  energy  consumption  reached  to  65.70  million  tons  of 
coal  equivalent.  Coal  consumption  accounts  for  30.3%  of  the  total 
energy  consumption.  According  to  the  statistic  data,  all  of  natural 
gas,  97%  of  coal,  80%  of  refined  oil,  and  70%  of  electricity  need  to 
be  transferred  from  the  neighboring  provinces  [61,62].  Particu¬ 
larly,  electricity  mainly  depends  on  “West-East  Electricity  Trans¬ 
mission  Project”,  and  an  implemented  project  named  “West-to- 
East  Gas  Transport”  helps  deliver  natural  gas  to  the  city.  The 
energy  redeployment  pattern  has  bonded  its  own  economic 
operation  and  the  external  supply  together,  which  could  arouse 
many  unpredictable  factors  and  increase  the  risk  of  economic 
development  in  Beijing.  The  city’s  power  demand  is  growing  with 
the  rapid  development  of  economic  and  continually  improvement 
of  people's  living  standard,  which  also  promotes  the  swift  devel¬ 
opment  of  it  power-generation  industry.  As  illustrated  in  Fig.  1, 
the  total  amount  of  electricity  consumption  in  Beijing  increased 
from  51.01  billion  kWh  in  2004  to  83.09  billion  kWh  in  2010. 
Elasticity  ratio  of  electricity  consumption  has  reached  more  than 
one  in  recent  years,  which  means  that  the  growth  rate  of 
electricity  demand  is  far  greater  than  the  speed  of  domestic 
economic  development.  Although  Beijing’s  energy  end-consump¬ 
tion  structure  has  improved  in  recent  years,  its  energy  consump¬ 
tion  is  still  relied  heavily  on  coal.  The  proportion  of  both  coal  and 
coke  occupies  40%,  and  75%  of  electricity  generated  from  coal- 
fired  conversion  technique.  The  proportion  of  coal  in  the  total 
energy  is  gradually  decreasing,  and  the  total  amount  of  coal  is 
continuing  to  increase.  Because  of  this,  a  large  amount  of  C02 
emission  was  caused  by  coal  consumption,  and  electric  power 


Fig.  1.  Annual  electricity  consumption  of  Beijing  from  2004  to  2009. 


plants  become  one  of  the  main  sources  of  C02  contribution  in 
China  [13,16]. 

Cun'ently,  the  city  still  faces  a  number  of  challenges  to  achieve 
the  C02-reduction  goals,  and  those  pressures  come  from  both 
domestic  and  overseas.  Firstly,  the  contradiction  between  electricity 
demands  and  C02-reduction  goals  would  be  increasingly  prominent. 
Secondly,  energy  structure  is  needed  to  further  adjust  and  optimize. 
Thirdly,  the  superiority  of  renewable  energy  has  not  been  fully 
played  by  now.  EPS  is  one  of  the  major  sectors  for  emission  C02. 
As  a  consequence,  dying  to  develop  carbon  trading  in  EPS,  stimulat¬ 
ing  the  enthusiasm  of  the  C02  reduction  through  market  means,  and 
reducing  system  cost  by  using  advanced  technology  will  vigorously 
promote  electric  power  industiy  restructuring  and  industrial  upgrad¬ 
ing.  Meanwhile,  China  has  been  viewed  to  be  one  of  the  most 
promising  emission  markets  for  reducing  GHG,  although  there  is  no 
emission  reduction  assignment.  This  study  attempts  to  introduce  CET 
into  the  Beijing’s  electric  power  market  to  deal  with  these  questions: 
(1)  how  to  determine  initial  distribution  of  permits  and  sector 
coverage  for  Beijing’s  domestic  power  market;  (2)  how  to  quantify 
the  final  system  cost  of  each  power  conversion  technology  under 
CET;  (3)  how  to  ascertain  the  amount  of  electricity  supply  and  C02 
mitigation.  Fig.  2  presents  the  CET  scheme  of  power  generation 
system  in  Beijing.  The  relationships  among  government,  fuel  market, 
power  producers,  carbon  market,  electricity  market,  economic  infor¬ 
mation  as  well  as  energy  flow  have  been  presented.  As  is  shown  in 
the  figure,  it  is  obvious  that  uncertainties  would  influence  the 
process  of  modeling  C02  emission  in  EPS  of  Beijing.  For  example, 
the  amount  of  C02  emission  generated  from  the  electricity  genera¬ 
tion  sector  may  vaiy  because  of  the  various  electricity  demands; 
distinctly,  for  such  a  large  system,  the  coefficients  related  to  cost  and 
benefit  information  is  not  sufficient.  An  integrated  energy  and 
environmental  management  system  can  be  characterized  by  one  or 
several  sources  (i.e.,  power  conversion  technology)  generally.  A  large 
amount  of  C02  emission  from  these  power  plants  may  lead  to 
adverse  impacts  on  climate  change.  For  example,  increasing  amount 
of  C02  in  the  atmosphere  may  affect  weather  condition  changes, 
sea/land  ice  cover  decreases,  biodiversity  changes,  and  ecosystem 
changes. 

3.2.  Modelling  formulation 

In  CET  management,  it  is  imperative  to  identify  desired  schemes 
for  energy-flow  allocation,  electricity-supply  and  C02  emission 
permits  with  minimized  system  cost  and  maximized  system  relia¬ 
bility.  Conventional  and  renewable  energy  resources  with  limited 
availabilities  are  employed  (including  raw  coal,  washed  coal,  coal 
products,  coal  oven  gas,  other  coal  gas,  diesel  oil,  fuel  oil,  refinery  dry 
gas,  natural  gas,  oil  products,  coking  products,  heating  power  and 
other  energy  sources).  Seven  power-conversion  technologies  are 
taken  into  consideration,  such  as  coal-fired  power,  hydropower, 
pumped  storage,  wind  power,  photovoltaic  power,  waste  generation 
and  biomass  power.  In  order  to  reduce  the  cost  for  C02  treatment, 
emission  trading  is  considered  for  each  power  conversion  technol¬ 
ogy.  A  target  quantity  of  C02  emission  quota  is  allocated  to  each 
power  conversion  technology.  If  this  quantity  is  satisfied,  the  EPS 
will  bring  net  benefits.  If  this  quantity  exceeds  the  regulated  level,  it 
is  necessary  to  take  measures  to  decrease  the  C02  emission.  In 
response  to  such  regulation,  the  power  plants  have  to  optimize  C02 
treated  to  achieve  a  minimized  system  cost  while  to  satisfy  the 
GHG-emission  requirement.  Through  the  CET  program,  each  power 
plant  can  sell  credit  to  other  power  plants  with  higher  electric 
power  profitability.  The  C02  emission  permits  can  thus  be  real¬ 
located  to  the  most  efficient  power  plants  instead  of  proportionally 
allocated  to  each  power  plant.  Consequently,  the  developed  FIMP 
method  can  be  used  to  address  these  problems,  leading  to  a  FIMP- 
based  carbon  emission  trading  (FIMP-CET)  model.  The  objective 
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Fig.  2.  Carbon  trading  in  electric  power  systems  of  Beijing. 


function  of  the  FIMP-CET  model  can  be  formulated  as  follows: 
Min/*  =  (1)  +  (2)  +  (3)  +  (4)  +  (5)  +  (6)  +  (7)  (8-1) 

(1)  Cost  for  purchasing  energy: 

E  E  CBNik  (VC )  x  ZHXik  X  XNL  *  (8  -  1  a) 

i  =  1  k  =  1 

(2)  Operating  cost  for  electricity  generation: 

EED»  ( PFDnk+ynk  X  A  PFD^+YPL^  x  YDAnkq)  x  ODC±(a±) 

n= 1 k= 1 q = 1 

(8 -lb) 


(3)  Penalty  for  electricity  shortage: 

E  Ph  x  WCC*  (v± )  X  OSE*  (8  -  lc) 

h  =  1 

(4)  Cost  for  capacity  expansion: 

E  E  E  L  X  YPLniq  x  YDAnkq  X  YD/±(y±  )  (8  -  Id) 

n=lfc=lq=l 

(5)  Cost  for  electricity  transmission: 

E  E  E  1  x  ( PFDnk+ynk  X  A PFDnk  +  YPLniq  X  YDAnkq )  x  CUT* 

n=\k=\q=\ 

(8-  le) 
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Table  1 

Purchase  cost  for  energy  source  (106  RMB  ¥/PJ). 


Time  period 

k= 1 

k= 2 

k= 3 

Raw  coal 

[22.90  (l+y,)',  25.80  (l+y,)'] 

[25.80  (l  +  y,)',  28.60  (l+y,)'] 

[28.90  (l+y,)',  32.50  (l  +  y,)'] 

Washed  coal 

[29.80  (l  +  y,)',  33.30(1+7,)'] 

[29.80  (l+y,)',  33.30  (l+y,)'] 

[29.80  (l+y,)',  33.30  (l+y,)'] 

Coal  products 

[29.80  (l+y,)',  29.80  (l+y,)'] 

[29.80  (l+y,)',  29.80  (l+y,)'] 

[29.80  (l+y,)',  29.80  (l+y,)'] 

Coal  oven  gas 

[29.80  (l+y,)',  29.80  (l+y,)'] 

[29.80  (1 +(»,)',  29.80  (l+y,)'] 

[29.80  (l+y,)',  29.80  (l+y,)'] 

Other  coal  gas 

[116.60  (l+y,)',  118.70  (l+y,)'] 

[118.60  (l+y,)',  122.60(1  +y,)'[ 

[119.80  (l+y,)',  125.70  (l+y,)'] 

Diesel  oil 

[143.30  (l+y,)',  148.00  (l+y,)'j 

[148.00  (l  +  y,)',  154.30  (l+y,)'] 

[156.20  (l+y,)',  138.50  (l+y,)'] 

Fuel  oil 

[131.10  (l+y,)',  133.50  (l+y,)'] 

[136.20  (l  +  y,)',  140.70  (l+y,)'] 

[142.10  (l+y,)',  84.40  (l+y, )') 

Refinery  dry  gas 

[29.80  (l+y,)',  29.80  (l+y,)'] 

[29.80  (l+y,)',  29.80  (l+y, )') 

[29.80  (l+y,)',  32.50  (l+y,)'] 

Natural  gas 

[46.40  (l+y,)',  50.00  (l+y,)'] 

[60.41  (l+y,)',  68.90  (l+y,)'] 

[74.40  (l+y,)',  32.50  (l+y,)'] 

Oil  products 

[123.40  (l+y,)',  128.20  (l+y, )'] 

[131.40  (l+y,)',  135.00  (l+y,)'] 

[135.00  (l+y,)',  138.50  (l+y,)'] 

Coking  products 

[36.80  (l+y,)',  50.00  (l+y,)'] 

[47.40  (l+y,)',  58.00  (l+y,)'] 

[56.10  (l+y,)'  84.40  (l+y,)'] 

Heating  power 

[30.00  (l+y,)',  47.40  (1  +y,)'[ 

[30.00  (1  +co,)',  32.00  (1  +y,)'] 

[30.00  (l+y,)',  32.50  (l+y,)'] 

Other  energy 

[64.14  (l  +  y,)',  32.00  (l+y,)'] 

[68.50  (1  +  co,)'.  71.90  (l+y,)'] 

[71.81  (l+y,)',  76.10  (l+y,)'] 

Note :  symbols  y,-  and  t  denote  the  interest  rate  and  time  interval  (from  the  time  of  data  sources  to  this  study). 


(6)  Cost  for  C02  reduction: 

E  E  E  L  x  ( PFDnk+ynk  X  &PFDnk  +  YPL±kq  X  YDAnkq ) 

n=lk=lq=l 

X (CPN±  +  SFk  X  CEE ±  -SU± )  (8  -  1  f ) 

(7)  Penalty  for  excess  C02  emission: 

E  E  E  L  x  Pi.  x  x  CMR  ±  (8  -  1  g) 

n  =  1  k  =  1  h  =  1 

where  / 1  is  total  cost  for  purchasing  energy,  electricity  generation, 
penalty  for  electricity  shortage,  capacity  expansion,  electricity 
transmission,  C02  reduction  and  penalty  for  excess  C02  emission;  i 
represents  type  of  purchasing  energy,  where  i=l  for  raw  coal;  i=  2 
for  washed  coal;  i=3  for  coal  products;  i=4  for  coal  oven  gas;  i=  5 
for  mixed  coal  gas;  i= 6  for  diesel  oil;  i= 7  for  fuel  oil;  /=  8  for 
refinery  dry  gas;  i= 9  for  natural  gas;  i=10  for  oil  products;  i=ll  for 
coking  product;  i=  12  for  heating  power;  i=  1 3  for  other  energy 
sources;  n  denotes  power  conversion  technology,  n=l  for  coal-fired 
power;  n= 2  for  hydropower;  n= 3  for  pumped  storage;  n= 4  for 
wind  power;  n=5  for  photovoltaic  power;  n  =  6  for  waste  genera¬ 
tion;  n=7  for  biomass  power;  k  is  planning  period,  k=  1,  2,  3;  L  is 
planning  time  period  (three  years);  q  is  capacity-expansion  pro¬ 
gram;  h  is  electricity-demand  level,  h=  1,  2,  3;  ph  is  probability  of 
demand-level  h  occurrence  in  period  k  (%);XNL^  is  amount  of 
energy  supply  from  source  i  in  period  k  (PJ);  OSE jy  is  the  amount 
of  electricity  shortage  (e.g.,  imported  electricity)  under  electricity- 
demand  level  h  (PJ);  XCP±k  is  the  amount  of  excess  C02  generated 
from  power  conversion  technology  n  in  period  k  (103  t);  YPL^k  is 
the  binary  variables  for  identifying  whether  or  not  a  capacity 
expansion  action  of  power  conversion  technology  n  under  different 
expansion  program  q  needs  to  be  undertaken  in  period  k;  PFD±k  is 
pre-regulated  electricity  target  of  via  power  conversion  technology  n 
which  is  promised  to  end-users  in  period  k  (PJ);  A  PFDnk  = 
PFD^k—PFD„k,  PFD^k  =  PFDnk+  A PFDnk  x  ynk,  ynk  e  [0,1];  CBNlk  is 
purchase  cost  for  energy  source  i  in  period  k  (106  RMB¥/PJ);  yjf  is 
interest  rate  in  period  k;  is  interest  rate  under  variable  electricity 
demand-level  h;  ZFIXik  is  transform  coefficient  (PJ/103  t  standard 
coal);  WCC^y  is  purchasing  cost  for  importing  power  from  other 
regions  under  variable  electricity  demand-level  h  (106  RMB¥/PJ); 
YDAqnk  is  vexpansion  capacity  for  conversion  technology  n  under 
expansion  program  q  in  period  k  (MW);  YDI~k  is  expansion  cost 
for  conversion  technology  n  in  period  k  (106  RMB¥/MW);  ODC±k 
is  operation  cost  for  power  conversion  technology  n  in  period 


k  (106  RMB¥/PJ);  is  electricity  equipment  depreciation  rate  for 
conversion  technology  n  in  period  k\  CUT f  is  cost  for  electricity 
transmission  in  period  k  (106  RMB¥/PJ);  CPNf  is  cost  for  reducing 
C02  emission  in  period  k  (106  RMB¥/PJ);  SFk  is  C02  reduction 
efficiency  in  period  /<;  CEEf  is  cost  for  C02  emission  in  period  k 
(106  RMB¥/PJ);  SU^  is  financial  subsidy  in  period  k  (106  RMB¥/PJ); 
GMR±k  is  operating  cost  for  excess  C02  released  from  power 
conversion  technology  n  in  period  k  (103  RMB¥/tonne). 

Table  1  reveals  the  purchase  cost  for  energy  sources.  In  energy 
systems,  energy  prices  are  closely  related  to  the  volatility  of  interest 
rates  and  time  intervals.  For  example,  for  long-term  EPS  planning, 
the  purchase  cost  for  energy  source  ( CBN )  can  be  affected  by  the 
fluctuation  of  interest  rates,  leading  to  that  the  lower  and  upper 
bounds  of  purchase  cost  can  vary  with  the  interest  rate.  Under  this 
situation,  the  concept  of  crisp  interval  may  not  be  suitable  for 
describing  this  uncertainty.  Functional  intervals  can  be  defined  as  a 
lower-  and  an  upper-bound,  which  are  both  functions  of  its 
associated  impact  factor,  can  be  used  to  effectively  reflect  such  an 
uncertainty.  For  example,  if  CBN^  is  expressed  as  a  functional 
interval  of  [22.90  (1  +  yi)t,  25.80  (1  +  yi)f]  x  106  RMB  ¥/PJ;  symbols  yi 
and  t  denote  the  interest  rate  and  time  interval  (from  the  time  of 
data  sources  to  this  study),  CBN^  is  a  function  of  interest  rate, 
ranging  between  22.90  (l+yi)'x  106  RMB¥/PJ  and  25.80 
(1  +"/'/)'  x  106  RMB¥/PJ.  Therefore,  the  definition  of  functional  inter¬ 
val  can  help  reflecting  modeling  uncertainties  with  more  complex¬ 
ities  way  and  describing  the  real-world  conditions  with  more 
efficiency  way.  The  systems  analysis  methods  can  be  able  to  reflect 
interactive,  complex,  dynamic  and  uncertain  features  of  the  CET 
management  systems.  The  outputs  would  be  interpreted  to  generate 
desired  planning  alternatives  for  a  number  of  human  activities,  as 
well  as  the  related  policies  and  strategies. 

The  system  cost  is  to  be  minimized  subject  to  a  set  of  constraints 
that  describe  various  impact  factors  and  their  interactions.  These 
constraints  contain  resource  availability  constraints,  mass  balance 
constraints,  power  demand-supply  balance,  C02  emission  con¬ 
straints,  reallocated  emission  permits  constraints,  expansion  capital 
constraints,  integer-variable  constraints,  and  nonnegative  con¬ 
straints.  These  constraints  include  a  number  of  inequalities  which 
define  relationships  among  various  decision  variables  and  system 
conditions,  which  can  be  formulated  as  follows: 

3.2.1.  Resource  availability  constraints 

13  3 

E  XNLit  ^  R20*  •  Yi.k  (8  -  2a) 

i  =  1  k  =  1 
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Table  2 

Electricity  demands  and  optimal  electricity  supply  scheme  (PJ). 


Probability  (%) 

Time  period 

k= 1 

k= 2 

k= 3 

Electricity  demands 

Low 

60 

[919.56,  5817.35] 

Medium 

20 

[1036.50,  6483.00] 

High 

Optimal  electricity  supply 

20 

[1145.22,  7081.32] 

Coal-fired  power 

[107.49,  111.79] 

[117.96,  122.68] 

[128.43,  133.56] 

Hydropower 

[30.14,  31.35] 

[30.00,  31.20] 

[38.52,  40.06] 

Pumped  storage 

[4.36,  4.54] 

[7.50,  7.81] 

[10.64,  11.07] 

Wind  power 

[7.54,  7.84] 

[8.79,  9.14] 

[10.05,  10.45] 

Photovoltaic  power 

[21.35,  22.21] 

[27.63,  28.74] 

[33.91,  35.27] 

Waste  generation 

[42.08,  43.76] 

[44.17,  45.94] 

[46.26,  48.111 

Biomass  power 

Solution  of  electricity  shortage 

[247.82,  257.74] 

[268.76,  279.51] 

[289.69,  301.28] 

Low 

60 

[207.31  479.22] 

Medium 

20 

[248.77,  525.01] 

High 

20 

[289.70,  579.81] 

where  RZO £  denotes  energy  availability  of  energy  type  i  in  period  fc 
(104 1  standard  coal).  This  constraint  represents  that,  for  each  power 
conversion  technology  in  each  period,  the  amount  of  energy  utili¬ 
zation  must  be  not  less  than  the  total  available  energy  amounts. 

3.2.2.  Mass  balance  constraints 


E  E  x  CZNi  *  E  E^PFDnk+ynk  x  Ap™nt),  vk.  n= 1 

i  =  1  k  =  1  ii  =  Ik  =  1 

(8 -2b) 

where  CZN £  is  the  transfer  coefficient  for  unit  coal  raw  con¬ 
sumption  for  electricity  production;  Constraint  (8-2b)  specifies 
that  the  amount  of  generated  electricity  in  power  conversion 
technology  must  not  exceed  its  existing  and  expanded  capacities. 
Mass-balance  constraint  is  established  to  ensure  that  the  input 
energy  is  greater  than  the  output  one. 

3.2.3.  Constraints  of  power  demand-supply  balance 


E  E  E  F  X  (PFF>nk  +Ynk  x  APFF)nk  +  YPFnkq  x  YF>Ankq 

n  =  lk=l<j  =  l 

xhk-SD±k  x  hk  +  OSEf )  x  (1  -SLk)  >  DEph,  Vh 

(8  -  2c) 

where  hk  is  electric  generation  hours  in  period  k  (hour);  SD±k  is 
the  power  lost  for  each  power  conversion  technology  in  each 
planning  period;  DEph  represents  random  variable  of  total  elec¬ 
tricity  demand  of  level  ft  in  period  k  (PJ).  Constraint  (8-2c) 
is  generated  for  each  sector  to  ensure  that  the  energy  outputs 
from  the  demand  technologies  be  equal  to  the  end-use  energy 
demands.  Table  2  depicts  the  electricity  demands  and  their 
associated  probabilities  of  occurrence.  According  to  Beijing  Sta¬ 
tistics  Bureau  (from  1991  to  2011),  three  discrete  target  values 
(i.e.,  low,  medium  and  high)  are  selected  as  the  range  of  interval. 
Additionally,  division  of  the  targets  into  a  number  of  predefined 
values  associated  with  probabilities  (20%,  60%,  and  20%)  can  also 
meet  the  requirement  of  two-stage  stochastic  programming  (TSP) 
[43,63]. 


3.2.4.  C02  emission  constraints 


EE  ELx  ( PPDnk+ynk  x  A  PFD-nk  +  YPLn\q  X  YDAnkq  x  hk) 

n=lk=lq=l 

xAMR  ±  x  (1  -SFk)  <  ESf  ,  Vk  (8  -  2d) 


7  3  3 

EEELx  ( PFDnk+ynk  X  APFDnk  +  YPLnkq  x  YDAnkq  X  fife) 

n=lk=lq=l 

xAMR±  x  (l-SFk)>XCP±k,  Vn.k  (8  —  2e) 

where  AMR±k  denotes  emission  amount  of  C02  for  power  conversion 
technology  n  in  period  k  (103  t/PJ);  ESf  is  the  allowed  amount  of  C02 
in  period  k  (103  t);  constraints  (8-2d)  and  (8-2e)  set  up  limitations  for 
the  total  amount  of  C02  emission  s;  they  denote  that  the  excess  C02 
generated  from  power  conversion  technology  should  conform  the 
allowed  amount  of  C02  formulated  by  the  government. 

3.2.5.  Reallocated  C02-emission  permits  constraints 


7  3  3 

EEELx  HPFDnk+ynk  X  APFDnk  +  YPL±q  X  YDAnkq  x  hk) 


n=lk=lq=l 

xAMR±k  -XCP„±]  < ZFX±  ,  Vn.k  (8  -  2f) 

EEE  ZFXnk  2:  (1  -ytp  )TFDni  .Vn.k  (8  -  2g) 

n=lk=lq=l 

ZFXn±  <  (1  -co)  x  TFD±  ,  Vn,k  (8  —  2h) 

ESf  =  Bf  x  l„  x  AMR  ±  x  DEf  ,  Vn.k  (8  -  2i) 


where  ZFX±k  is  the  reallocated  C02  emission  permit  generated 
from  power  conversion  technology  n  in  period  k  (103t);  TFD±k 
represents  the  discharge  limit  of  total  C02  emissions  for  power 
conversion  technology  n  in  period  k  (103t);  co  is  C02-reduction 
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efficiency  level  (percentage  of  reduced  total  C02  emission  per¬ 
mit);  Bf  is  power  demand  and  supply  index  in  period  fc;  l„  is 
amount  of  C02  emission  loading  per  PJ  electricity  for  power 
conversion  technology  n.  In  Beijing,  each  power  conversion 
technology  has  a  pre-regulated  generation  target.  If  the  target  is 
not  exceeded,  the  system  will  be  encountered  the  regular  cost; 
otherwise,  the  system  will  be  subject  to  penalties  resulted  from 
the  extra  labor,  management,  operation  and  maintenance  costs, 
or  capacity  expansion  and  higher  costs  for  importing  clean 
energy.  For  example,  in  power  plant,  the  amounts  of  C02  emis¬ 
sions  vary  qualitatively  and  quantitatively  from  one  power  plant 
to  another,  which  can  result  in  huge  variations  in  the  cost  of 
achieving  targets  of  emission  limits.  This  difference  in  cost  can 
also  encourage  managers  of  power  plants  to  carry  out  carbon 
emissions  trading  scheme.  In  order  to  reduce  the  cost  for  C02 
treatment,  CET  is  considered  for  all  the  power  plants.  Based  on  the 
local  C02  emission  management  policies,  a  target  quantity  of  C02 
emission  quota  is  allocated  to  each  power  plant.  If  this  quantity  is 
satisfied,  the  power  generation  system  will  bring  net  benefit.  If 
this  quantity  exceeds  the  regulated  level,  power  plants  will  have 
to  take  measures  to  decrease  the  C02  emission.  In  response  to 
such  regulation,  the  power  plants  need  to  optimize  C02  treated  to 
achieve  a  maximized  system  net  benefit  while  to  satisfy  the  GHG 
emission  requirement.  Through  the  CET,  each  power  plant  can  sell 
credit  to  other  power  plants  with  higher  electric  power  profit¬ 
ability.  The  C02  emission  permits  can  thus  be  reallocated  to  the 
most  efficient  power  plants  instead  of  proportionally  allocated  to 
each  power  plant.  From  a  long-term  planning  view,  decision¬ 
makers  have  to  face  a  dilemma  of  either  investing  more  funds  on 
capacity  expansion  of  existing  facilities  or  turning  to  other  energy 
production  options  or  putting  extra  funds  into  energy  imports  at 
raised  prices  [50].  These  constraints  about  reallocated  C02-emis- 
sion  permits  are  used  for  guarantee  the  implement  of  CET. 

3.2.6.  Expansion  capital  constraints 

HiLx  YPL  niq  x  YDAnkq  X  YD/±(-4)  <  MN±  Vn.fc 

n=lfc=lq=l 

(8  -  2j) 

where  MN±k  denotes  the  expense  for  power  plant  expansion 
for  power  conversion  technology  n  in  period  k  (106  RMB¥);  Constraint 
(8-2j)  is  formulated  to  secure  sufficient  financial  capacities  for 
satisfying  production  of  a  given  energy,  where  expansion  constraints 
describe  the  lower  and  upper  bounds  of  capacity  expandability. 

3.2.7.  Interval-integer  variable  constraints 


YPL 


nkq 


{ 


=  1;  if  capacity  expansion  is  undertaken 
=  0;  if  otherwise 


Vn.fc,  q 


(8  -  21<) 


From  a  long-term  planning  point  of  view,  the  city’s  electricity 
demands  keep  increasing  due  to  population  increase  and  eco¬ 
nomic  development.  This  tendency  can  lead  to  insufficient  capa¬ 
city  to  meet  the  overall  electricity  demands  from  multiple  end- 
users.  The  related  optimization  analysis  will  require  the  use  of 
integer  variables  to  indicate  whether  a  particular  facility  devel¬ 
opment  or  expansion  option  needs  to  be  undertaken  [63]. 


3.2.8.  Nonnegative  constraints  for  decision  variables 


XNLf  ,OSE±  ,XCP±,ZFX±  >  0 


(8-21) 


The  above  model  can  be  solved  according  to  the  solution 
method  as  described  in  Section  2.  In  the  above  modeling  for¬ 
mulation,  pre-regulated  electricity  target  (PFD^k )  are  expressed  as 
interval  numbers.  An  optimized  set  of  PFD±k  values  will  be 
identified  by  having  ynk  as  decision  variables;  this  optimized  set 
may  correspond  to  a  minimized  system  cost.  In  detail,  let 
PFDnk  =PFDn, k+Ynk  x  APFD±  ,  where  A PFDnk  =  PFDfk-PFDnk  and 
ynk  e  [0,1].  As  a  result,  when  PFD^  approach  their  lower  bounds 
(i.e„  when  ynk  =  0),  a  relatively  low  cost  would  be  obtained; 
however,  a  higher  penalty  may  have  to  be  paid  when  the 
electricity  demand  is  not  satisfied.  Conversely,  when  PFD^  reach 
their  upper  bounds  (i.e.,  when  ynk=  1),  a  higher  cost  would  be 
generated  but,  at  the  same  time,  a  lower  risk  of  violating  the 
promised  targets.  The  objective  function  [i.e.,  Eq.  (8-1)]  is  to 
minimize  the  expected  system  cost.  If  the  promised  electricity 
target  is  delivered,  a  net  cost  to  the  local  economy  will  be 
generated  for  each  unit  of  electricity  supply.  However,  if  the 
promised  electricity  is  not  delivered,  either  the  electricity  must 
be  obtained  from  higher-priced  alternatives,  or  the  demands  must 
be  curtailed  with  the  costs  of  reduced  industrial  and/or  agricul¬ 
tural  productions.  These  would  then  result  in  an  increase  of 
system  cost  to  each  user.  Fig.  3  presents  the  framework  of 
FIMP-CET  model  of  Beijing’s  EPS.  In  this  system,  different  uncer¬ 
tainties  are  suffused  in  input  data  including  price  data,  demand 
data,  emission  data  and  availability  data.  By  solving  the  model, 
the  schemes  about  energy  supply,  electricity  supply,  carbon 
quotas  and  capacity  expansion  would  be  generated.  Moreover, 
excess  and  reallocated  C02  can  be  achieved  when  the  allowable 
C02  levels  (as  pre-regulated  by  authorities)  exceed. 


4.  Result  analysis 

In  this  study,  three  time  periods  are  considered,  with  each 
having  an  interval  of  3  years.  Over  the  9-year  planning  horizon,  an 
existing  power  generation  system  is  available  to  meet  electricity 
needs  of  the  city.  Facilities  of  coal-fired  power  plant,  hydropower 
plant,  pumped  storage  plant,  wind  power  plant,  photovoltaic 
power  plant,  waste  generation  plant  and  biomass  power  plant 
are  available  for  power  generation.  The  results  may  shed  light  on 
the  advantages  of  the  government’s  decision  and  even  potentially 
cause  a  change  in  policy. 

4.1.  Carbon  trading  analysis 

In  this  study,  CET  is  considered  for  helping  satisfy  the  GHG 
emission  requirement  of  EPS.  Based  on  local  C02  emission 
management  policies,  trading  amount  of  C02-emission  permit 
would  be  allocated  to  each  power  plant.  If  C02  emission  equal  to 
(or  less  than)  the  pre-regulated  quota,  benefits  would  be  achieved. 
In  comparison,  if  C02  emission  exceeds  the  pre-regulated  quota, 
reduction  measures  should  be  adopted  in  power  plants  to 
decrease  the  C02  emission  for  avoiding  high  system  cost.  To 
response  to  such  regulation,  the  power  plants  need  to  generate 
trading  amount  of  C02-emission  permit  to  achieve  a  minimized 
system  cost  while  to  satisfy  the  GHG  emission  requirement.  Fig.  4 
depicts  the  solutions  of  amount  of  CET  by  solving  the  FIMP-CET 
model.  Coal-fired  power  plant  would  purchase  partial  C02-emis- 
sion  permit  from  other  power  plants  (e.g.,  hydropower,  pumped 
storage,  wind  power,  photovoltaic  power,  waste  generation  and 
biomass  power)  over  the  planning  horizon.  For  example,  the 
amounts  of  C02  emissions  purchased  from  other  power  plants 
would  be  [4186.63,  5116.12]  x  103 1  in  period  1,  [4159.19, 
5067.86]  x  103  t  in  period  2,  and  [4145.48,  5019.59]  x  103  t  in 
period  3.  Sellers  would  be  the  roles  for  hydropower,  pumped 
storage,  wind  power,  photovoltaic  power,  waste  generation  and 


C02trading  quantity  (103t) 
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Fig.  3.  Framework  of  FIMP-CET  model  of  Beijing. 


■  Coal-fired  power  ■  Hydropower  ■  Pumped  storage  ■  Wind  power 

■  Photovoltaic  power  □  Waste  generation  ■  Biomass  power 

Fig.  4.  Solution  of  trading  amount  of  C02-emission  permit. 


biomass  power  plants.  For  example,  for  the  hydropower  plant, 
the  amount  of  C02-emission  permit  sold  to  other  power  plants 
(e.g.,  coal-fired)  would  be  [1930.43,  2152.26]  x  103  t,  [1930.43, 
2131.95]  xl03t  and  [1930.43,  2111.65]  x  103  t  during  periods 
1  to  3.  For  the  biomass  power  plant,  the  amounts  of  CET  would 
be  [816.56,  971.51]  x  103,  [816.56,  962.34]  x  103  and  [816.56, 
953.18]  x  103  t  during  periods  1  to  3.  For  the  wind  power  plant, 
the  amounts  would  be  [134.42,  290.79]  xl03t  in  period  1, 
[106.99,  288.04]  x  103  t  in  period  2,  and  [93.27,  285.30]  x  103  t 
in  period  3.  So,  C02  emission  permits  can  be  reallocated  to  the 
most  efficient  power  plants  instead  of  proportionally  allocated  to 
each  power  plant. 
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4.2.  Energy  resources  supply 

Energy  resources  used  in  Beijing  can  be  further  divided  into 
seven  types,  such  as  raw  coal,  washed  coal,  coal  products,  coal 
oven  gas,  other  coal  gas,  diesel  oil,  fuel  oil,  refinery  dry  gas, 
natural  gas,  oil  products,  coking  products,  heating  power  and 
others.  Table  3  reveals  the  results  of  optimal  energy  supply 
(XNL;J  )  for  ensuring  the  normal  operation  of  EPS.  According  to 
the  table,  the  amount  of  raw  coal  would  decrease  from  period  1  to 
period  3.  For  example,  the  amount  of  raw  coal  would  be  [1668.24, 
1868.43]  PJ  in  period  1,  [1526.54,  1724.99]  PJ  in  period  2,  and 
[1426.54,  1640.52]  PJ  in  period  3.  Such  a  drop  is  due  to  closure 
small-scale  mining  facilities  for  reaching  stringent  environmental 
requirement.  The  gap  would  be  filled  by  imports  from  other 
energy  systems,  mostly  from  the  Shanxi  and  Inner-Mongolia 
provincial  energy  systems.  Meanwhile,  for  declining  high  external 
dependence  of  the  energy  utilization,  the  amount  of  diesel  oil  and 
fuel  oil  would  distinctly  decrease  from  period  1  to  period  3. 
Natural  gas  has  a  dramatic  growth  among  the  whole  planning 
horizon.  For  instance,  the  amount  of  natural  gas  would  reach 
[393.66,  440.90]  PJ  in  period  1,  [501.74,  566.97]  PJ  in  period  2,  and 
[602.44,  692.81]  PJ  in  period  3.  Thus,  natural  gas  would  be  one  of 
the  major  energy  sources  in  the  future  for  the  city.  Consequently, 
the  city’s  energy  policies  for  enhancing  energy  supply  capacity, 
speeding  up  the  utilization  of  renewable  energy,  and  increasing 
energy  infrastructure  investment  are  desired. 

4.3.  Electricity  supply 

Solutions  of  optimal  electricity  supply  scheme  can  be  inter¬ 
preted  based  on  the  results  presented  in  Table  2.  Several  power 
conversion  technologies,  including  coal-fired  power,  hydropower, 
pumped  storage,  wind  power,  photovoltaic  power,  waste  genera¬ 
tion  and  biomass  power.  Managers  would  have  to  promise  pre¬ 
regulated  target  values  of  electricity  supply  in  order  to  satisfy  the 
electricity  demand  objective.  Among  the  power  conversion  tech¬ 
nology,  optimized  electricity-generation  target  can  be  obtained. 
For  example,  biomass  power  would  become  the  most  competitive 
form  of  power  generation  in  the  near  future.  The  amount  of 
biomass  power  would  increase  from  [247.82,  257.74]  PJ  in  period 
1  to  [289.69,  301.28]  PJ  in  period  3.  The  raise  of  biomass  power 
generation  would  occur  in  whole  planning  horizon,  which  has  an 
increase  from  8.44%  to  8.45%.  For  environmental  consideration, 
the  generation  amount  of  coal-fired  power  will  be  tightened  in  a 
long-term  period.  For  example,  the  generation  amount  of  coal- 
fired  power  would  be  [107.49,  111.79]  PJ  in  period  1,  [117.96, 
122.68]  PJ  in  period  2,  and  [128.43,  133.56]  PJ  in  period  3, 


Table  3 

Solution  of  optimal  energy  supply  (PJ). 


Time  period 

k= 1 

k= 2 

k= 3 

Raw  coal 

[1668.24,  1868.43] 

[1526.54,  1724.99] 

11426.54,  1640.52] 

Washed  coal 

[7.58,  8.48] 

[9.23,  10.42] 

[11.23,  12.91] 

Coal  products 

[10.75,  12.04] 

[7.20,  8.14] 

[7.10,  8.17] 

Coal  oven  gas 

[1.93,  2.17] 

[2.29,  2.58] 

[2.49,  2.86] 

Other  coal  gas 

[34.94,  39.14] 

[34.37,  38.84] 

[34.35,  39.50] 

Diesel  oil 

[0.65,  0.73] 

[0.44,  0.49] 

[0.45,  0.52] 

Fuel  oil 

[10.41,  11.66] 

[3.43,  3.87] 

[4.22,  4.85] 

Refinery  dry  gas 

[2.07,  2.32] 

[3.91,  4.42] 

|3.75,  4.31] 

Natural  gas 

[393.66,  440.90] 

[501.74,  566.97] 

1602.44,  692.81] 

Oil  products 

[5.67,  6.35] 

[6.10,  6.89] 

|6.10,  7.02] 

Coking  products 

[27.39,  30.68] 

[22.90,  25.87] 

[20.33,  23.38] 

Heating  power 

[34.67,  38.83] 

[58.09,  65.64] 

[67.80,  77.97] 

Other  energy 

[14.70,  16.46] 

[15.60,  17.63] 

[19.00,  21.85] 

respectively.  What’s  more,  by  contrast  with  wind  power,  photo¬ 
voltaic  power  and  waste  generation  would  have  relatively  wide 
development  space.  The  amount  of  photovoltaic  power  grew 
steadily  in  the  planning  period,  changing  from  [21.35,  22.21]  PJ 
in  period  1  to  [33.91,  35.27]  PJ  in  period  3.  Moreover,  the  amount 
of  electricity  generated  from  waste  power  plant  would  go  up  with 
a  dramatic  rise  from  [42.08,  43.76]  PJ  in  period  1  to  [46.26, 
48.11]  PJ  in  period  3.  This  is  a  sharp  increase  of  9.93-9.94%.  In 
addition,  with  the  rapid  economic  growth,  electricity  demand  is 
increasing  year  by  year.  Electricity  generated  from  local  power 
plant  cannot  meet  such  a  huge  demand  for  electricity  consump¬ 
tion,  and  then  electricity  shortage  (O SEf  )  would  occur.  Thus,  for 
filling  the  electricity  demand  gap,  electricity  would  be  imported 
from  other  regions,  mostly  from  the  Shanxi  and  Inner-Mongolia 
provincial  energy  systems.  High  electricity  demand  and  the  total 
amount  of  C02  emitted  would  be  confined  with  a  certain  level 
during  the  planning  periods.  The  amounts  of  imported  electricity 
would  be  [207.31,  217.67]  PJ  in  period  1,  [248.77,  261.21]  PJ  in 
period  2,  and  [289.69,  304.18]  PJ  in  period  3.  The  imported 
electricity  occupied  7.61-18.40%  of  the  total  electricity  generation 
under  low  electricity-demand  level,  7.49-19.36%  under  medium 
electricity-demand  level,  and  7.57-20.19%  under  high  electricity- 
demand  level.  Although  electricity  demand  is  huge,  renewable 
energy  generation  contains  enormous  development  potential. 
It  is  requisite  to  promote  the  utilization  of  renewable  energy. 
As  important  measures  to  guarantee  a  nation’s  energy  security, 
the  development  and  utilization  of  renewable  energy  are  crucial 
for  supporting  C02  emission  reduction. 

4.4.  Capacity  expansion 

In  EPS,  existing  capacity  always  cannot  meet  the  electricity 
demands;  as  a  result,  facility  expansion  would  be  obtained  through 
the  F1MP-CET  model  to  guarantee  the  electricity  requirement  use. 
Table  4  expresses  the  results  of  expansion  schemes  (YPL±nk)  and 
capacity  investment  {MN±k ).  Capacity  expansion  in  this  study  means 
to  expand  production  scale  of  the  existing  power  plant,  to  eliminate 
backward  high  energy-consumption  equipment,  and  to  promote 
energy-saving  technologies  and  equipment.  In  period  1,  hydropower 
capacity  would  be  expanded  to  [80,90]  MW;  an  additional  capacity  of 
[80,90]  MW  from  wind  power  facilities  would  be  installed,  and  the 
amount  of  capacity  expansion  for  biomass  power  would  reach  to 
[80,90]  MW  as  well.  In  period  2,  hydropower  capacity  would  be 
expanded  to  [120,130]  MW,  and  an  additional  capacity  of  [120,130] 
MW  from  wind  power  facilities  would  achieve  installed.  In  period  3, 


Table  4 

Solution  of  expansion  scheme  and  capacity  investment. 


Time  period 

k=  1 

k=2 

k= 3 

Expansion  schemes 

(MW) 

Coal-fired  power 

[0,0] 

[0,0] 

[0,  0] 

Hydropower 

[80,90] 

[120,130] 

[0,  0] 

Pumped  storage 

[0,  0] 

[0,  0] 

10,  0] 

Wind  power 

[80,90] 

[120,130] 

[0,  0] 

Photovoltaic  power 

10, 0] 

[0,  0] 

[150,160] 

Waste  generation 

[0,0] 

[0,0] 

[0,  0] 

Biomass  power 

[80,90] 

[0,0] 

[120,130] 

Capacity  investment  (109  RMB¥) 

Coal-fired  power 

[0,0] 

[0,0] 

[0,  0] 

Hydropower 

[11.76,  25.83] 

[21.89,  37.32] 

10,  0] 

Pumped  storage 

10,  0] 

[0,  0] 

[0,  0] 

Wind  power 

[104.52,  246.04] 

[165.89,  284.31] 

[0,  0] 

Photovoltaic  power 

[0,0] 

[0,0] 

[1916.01,  1749.60] 

Waste  generation 

[0,0] 

[0,0] 

[0,  0] 

Biomass  power 

[53.57,  110.72] 

[0,0] 

[139.35,  140.38] 
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the  amount  of  capacity  expansion  for  photovoltaic  power  would  be 
[150,160]  MW,  and  an  additional  capacity  of  [120,130]  MW  from 
biomass  power  facilities  would  achieve  installed  as  well.  In  addition, 
capacity  investment  expenses  related  to  interest  rate  were  generated 
from  the  proposed  model  as  shown  in  table.  As  a  matter  of  fact, 
backward  equipment  would  cause  material  waste,  environmental 
pollution,  as  well  as  GHG  emissions.  As  a  result,  after  renewable 
energy  generation  expansions  occur  in  the  planning  periods,  the  total 
capacity  of  the  facilities  would  be  sufficient  to  meet  the  increased 
electricity  demand  and  the  reduced  C02  emission  requirement. 

4.5.  System  cost 

Under  multiple  uncertainties  in  the  related  factors  and  para¬ 
meters,  solution  of  the  objective  function  values  (/0pt  =  RMB¥ 
[173.51,  237.49]  x  1012)  represents  two  extremes  of  net  system  cost 
over  the  planning  horizon.  Conventionally,  as  the  actual  value  of  each 
continuous  variable  varies  within  its  lower  and  upper  bounds,  the  net 
system  cost  would  change  correspondingly  between  f~pt  and  fp  t. 
Normally,  a  plan  with  a  lower  system  cost  would  correspond  to  a 
lower  end-user  demand  and  a  lower  level  of  projection  to  economic 
development  and  population  growth  in  the  study  region;  however,  it 
would  result  in  a  higher  risk  of  violating  the  system  constraints. 
Besides,  the  upper  bounds  of  cost  coefficients  correspond  to  fp  t, 
implying  that  the  decision  maker  has  an  extremely  conservative 
attitude  for  estimating  the  system  cost  and  end-user  demand,  leading 
to  a  low  system-failure  risk  and  high  reliability  level  for  energy 
production  and  supply.  Conversely,  a  plan  with  a  lower  system  cost 
would  correspond  to  an  optimistic  estimation  towards  energy 
demand  projection  and  relevant  costs  (i.e„  lower-bound  system  cost, 
lower  energy  demand  and  lower  cost  coefficients),  which  also  means 
that  conservative  prediction  of  economic  development  and  popula¬ 
tion  growth  over  the  planning  horizon.  Hence,  a  decision  with  lower 
system  cost  could  lead  to  a  higher  risk  of  system  failure  and  lower 
system  reliability  for  meeting  energy  needs  in  the  study  region.  This 
demonstrates  that  there  is  a  tradeoff  between  system  cost  and 
reliability  level  for  energy  resources  management.  Table  5  represents 
the  solutions  for  system  cost,  which  is  composed  of  expenses  for 
energy  resources  purchase,  importing  power,  electricity  generation, 
capacity  expansion,  C02-emission  reduction,  and  penalty  for  excess 
C02  emission.  In  practice,  planning  for  the  high-bound  system  cost 
would  lead  to  a  lower  risk  of  violating  the  allowable  C02  emission 
level.  Conversely,  planning  for  the  low-bound  system  cost  would  lead 
to  a  higher  probability  of  violating  the  allowance.  Therefore, 
there  is  a  tradeoff  between  the  system  cost  and  C02  emission- 
allowance  violation  risk. 

4.6.  C02-reduction  efficiency 

CET  is  one  of  the  measures  to  achieve  the  C02  emission 
reduction  targets.  However,  it  is  essential  to  decrease  C02  emis¬ 
sion  by  using  advanced  technologies  in  each  power  plant.  Accord¬ 
ingly,  it  is  critical  to  discuss  the  complicated  relationship  between 


system  cost  and  various  C02-reduction  efficiency  level  ( SFk )  in 
response  to  CET  planning  in  EPS  of  Beijing.  By  solving  F1MP  model, 
it  is  obvious  that  the  optimized  system  costs  would  have  changes 
with  the  fluctuation  of  C02-reduction  efficiency  level.  Fig.  5 
declares  the  optimized  system  costs  of  the  model  under  different 
C02-reduction  efficiency  level.  There  would  be  a  downward 
tendency  of  optimized  system  cost  over  the  planning  horizon. 
For  example,  the  system  cost  would  decrease  from  RMB¥  [309.28, 
385.01]  xlO12  when  SFk  =  0%  to  RMB¥  [211.07,  275.55]  x  1012 
when  SFk  =  60%.  Then,  the  graph  throws  a  new  light  on  the 
situation.  For  example,  the  system  cost  would  increase  from 
RMB¥  [215.36,  279.77]  x  1012  when  SFk  =  70%  to  RMB¥  [223.93, 
288.53]  x  1012  when  SFk= 90%.  It  is  obviously  that  C02  emission 
permits  for  each  power  plant  would  be  optimized  through  C02 
emission  permits  trading  scheme,  and  the  trading  scheme  would 
be  more  effectual  under  the  condition  that  the  C02-reduction 
efficiency  level  equals  to  60%. 

Excess  C02  treated  by  different  power  plants  would  be  various 
under  different  C02-reduction  efficiency  levels,  as  shown  in  Fig.  6. 
For  example,  the  amount  of  excess  C02  treated  at  cold-fired 
power  plant  would  decrease  from  [29.31,  30.60]  xl06t  when 
SFk= 0%  to  [15.10,  15.58]  x  106t  when  SFk= 60%.  From  reduction 
efficiency  level  reaching  to  60%,  there  was  a  gradual  tendency  for 
excess  C02  treated  at  coal-fired  power  plant.  The  amount  of 
excess  C02  treated  at  coal-fired  power  plant  would  be  the  same 
([15.10,  15.58]  xl06t)  when  C02-reduction  efficiency  level 
changes  from  60%  to  90%.  When  C02-reduction  efficiency  level 
equals  to  60%,  the  excess  C02  treated  in  cold-fired  power  plant 
would  be  achieved.  Otherwise,  Fig.  7  describes  the  excess  C02 
treated  at  hydropower  plant.  Comparing  to  coal-fired  power 
plant,  hydropower  plant  would  have  a  low  level  of  excess  C02 
treated.  As  the  reduction  efficiency  increasing,  there  is  a  steadily 
decreasing  tendency  of  the  amount  of  excess  C02  treated  at 
hydropower  plant.  For  example,  the  amount  of  excess  C02  treated 
in  hydropower  plant  would  decline  from  [27.15,  59.02]  xl03t 
when  SFk= 0%  to  [2.71,  5.90]  x  103  t  when  SFk  =  90%.  It  indicated 
that  the  higher  the  reduction  efficiency  level  reached,  the  lower 


Fig.  5.  System  cost  under  different  reduction  efficiency  level. 


Table  5 

Each  activity  expense  1 09  RMB¥. 


Time  period 

k= 1 

k= 2 

k= 3 

Cost  for  purchasing  energy 

Cost  for  capacity  expansion 

Operating  cost  for  electricity  generation 

Cost  for  electricity  transmission 

Cost  for  C02  reduction 

Penalty  for  excess  C02  emission 

[2787.07,  5736.38] 

[729.06,  1476.23] 

[584.31,  630.63] 

[1.16,  1.25] 

[8404.07,  9149.43] 

[22926.42,  24252.99] 

[4121.44,  14307.79] 

[1087.49,  3493.00] 

[580.83,  635.43] 

[1.21,  1.31] 

[9389.17,  10211.28] 

[82163.52,  87565.08) 

[9047.31,  21734.193] 
[634.52,  1394.01] 
[552.72,  633.02] 

[1.28,  1.38] 

[10581.28,  11497.41] 
[41038.76,  44148.54] 
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Fig.  6.  Excess  C02  treated  in  cold-fired  power  plant. 


Fig.  7.  Excess  C02  treated  in  hydropower  plant,  (a)  Upper  bound  and  (b)  Lower 
bound. 


the  amount  of  excess  C02  which  should  be  treated.  This  is  due  to 
the  following  facts:  when  the  reduction  efficiency  level  is  low,  the 
C02  emission  permits  allocated  to  each  power  plant  are  relative 
low.  So,  less  C02  surplus  from  each  power  plant  would  be  treated, 
and  low  system  cost  would  be  gained.  However,  these  scenarios 
may  lead  to  lower  C02  reduction  level. 


5.  Discussion 

The  F1MP-CET  model  provides  a  linkage  to  pre-regulated 
management  policies  that  have  to  be  followed  when  a  modeling 
effort  is  undertaken.  If  not  considering  CET  planning,  it  can  be 
converted  into  a  FIMP-based  electric  power  systems  planning 
(FIMP-EPS)  problem.  The  solution  of  /±t  =  RMB¥  [196.13, 
262.08]  x  1012  from  FIMP-EPS  indicates  that  the  obtained  system 
cost  is  higher  than  that  from  the  FIMP-CET  model.  The  system 
cost  under  trading  scheme  is  obviously  lower  than  that  under 
non-trading  scheme.  Since  the  emission  trading  factors  are  not 
addressed  in  the  FIMP-EPS  model,  the  relevant  constraints  are 


relaxed.  Moreover,  in  order  to  meet  the  discharge  limit  of  total 
C02  emissions,  excess  C02  should  be  treated  by  other  measures, 
such  as  capture  and  storage,  chemical  absorption  etc.,  which  may 
further  increase  the  system  cost.  In  comparison,  considering  CET, 
the  FIMP-CET  model  can  effectively  reduce  system  cost.  To  ensure 
the  constraints  to  be  satisfied  under  all  interest  rate  levels, 
stringent  constraints  are  used  in  FIMP-CET,  although  this  might 
raise  the  system  cost.  In  FIMP-CET,  C02-emission  permits  for  each 
power  plant  are  optimized  through  C02  emission  permits  trading 
scheme,  and  the  trading  scheme  would  be  more  effectual  under 
the  condition  that  the  C02-reduction  efficiency  is  around 
60%.  This  may  also  lead  to  the  changing  of  the  system  cost  for 
EPS.  Therefore,  CET  is  effective  for  C02  permit  reallocation, 
and  different  policies  for  C02  management  are  associated  with 
different  levels  of  C02  management  cost  and  C02  mitigation- 
failure  risk. 

Compared  with  the  previous  studies,  the  FIMP-CET  model  has 
the  following  advantages.  Firstly,  the  FIMP-CET  has  advantages  in 
uncertainty  reflection  and  policy  analysis,  particularly  when  the 
input  parameters  are  provided  as  crisp  intervals,  functional 
intervals  and  probabilistic  distributions.  The  conventional  IPP 
method  can  only  tackle  intervals  with  fixed  lower  and  upper 
bounds.  In  the  FIMP,  constraints  can  be  adjusted  in  response  to 
the  fluctuations  of  external  factors.  The  infinite  objectives  can  be 
converted  into  a  single  one  and  the  infinite  constraints  can  also  be 
transformed  into  finite  ones.  Secondly,  this  is  the  first  attempt  to 
introduce  CET  into  Beijing’s  EPS  to  mitigate  C02  emissions 
through  the  developed  FIMP  method.  Electric  power  industry  is 
one  of  the  major  sources  of  C02  emission  in  the  City  of  Beijing. 
It  is  necessary  to  accumulate  relevant  experience  to  provide  a 
reliable  basis  for  establishing  a  regional  CET  market..  The  model 
also  incorporates  existing  international  CET  policies  directly  into 
optimization  process.  It  improves  upon  the  existing  approaches 
for  managing  C02  emissions  with  trading  scheme  of  EPS,  such 
that  robustness  of  the  optimization  process  can  be  enhanced. 
Third,  through  managing  C02  emissions  with  trading  scheme  of 
EPS  in  Beijing,  cost-effective  options  are  obtained  based  on  a 
least-cost  strategy  with  a  minimized  economic  cost  and  a  max¬ 
imized  electricity-supply  security.  The  obtained  solutions  provide 
more  practical  decision  bases  for  formulating  C02-reduction 
policies  and  assessing  the  associated  economic  implications 
in  purchasing  emission  permits  or  bearing  economic  penalties. 
The  continuous  variable  solutions  are  related  to  decisions  of 
energy  supply,  electricity  supply,  carbon  quota  and  capacity 
expansion,  the  interval  solutions  can  help  decision  makers  obtain 
multiple  decision  alternatives,  and  the  binary-variable  results 
represent  the  decisions  of  facility  expansion.  Thus,  decision 
makers  can  adjust  the  existing  demand  and  supply  patterns  of 
energy  resources,  generate  facilitate  dynamic  analysis  for  capacity 
expansion,  and  coordinate  the  conflict  interactions  among  eco¬ 
nomic  cost,  system  efficiency,  pollutant  mitigation  and  electricity- 
supply  security,  which  can  help  analyzing  various  policy  that 
are  related  to  different  levels  of  economic  penalties  when  the  pre¬ 
regulated  electricity  generation  schemes  are  violated. 

There  are  several  assumptions  for  formulating  the  FIMP-CET 
model,  which  may  bring  some  limitations  for  managing  C02 
emissions  with  trading  scheme  of  EPS  in  Beijing.  First,  normal 
distribution  is  assumed  to  the  availability  of  electricity  demands 
for  the  city;  in  probability  theory,  according  to  law  of  large 
numbers  and  central  limit  theorem,  any  random  event  approx¬ 
imates  a  normal  distribution  when  the  samples’  number  is  greater 
than  a  certain  number.  Second,  capacity  expansion  of  each  power 
conversion  technology  is  limited  to  once  within  the  planning 
horizon;  this  is  due  to  the  fact  that  the  expansion  capacity  of 
each  power  conversion  technology  is  restricted  by  the  condition 
of  finance  investment  and  facility  service  life.  Third,  each 
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community  in  the  city  has  the  same  economy  and  energy 
structures;  if  each  community  is  considered  as  different  economy 
and  energy  structures  such  that  more  complexities  would  be 
generated  which  may  beyond  the  model’s  mathematical  expres¬ 
sion  capacity. 

Besides,  techniques  of  post-optimality  analysis  (e.g.,  multi¬ 
criteria  decision  analysis,  analytical  hierarchy  process  technique, 
dual  programming,  and  parametric  programming)  could  be  used 
for  further  supporting  fine  adjustments  of  the  modeling  results 
and  thus  for  enhancing  their  applicability  to  practical  situations 
[64, 65]. Furthermore,  intelligent  decision  support  system  (IDSS) 
could  be  developed  based  on  an  integration  of  optimization 
modeling,  scenario  development,  user  interaction,  policy  analysis 
and  visual  display  into  a  general  framework.  Uncertainties  in  CET 
systems  could  be  effectively  reflected  and  addressed  through  the 
full-infinite  interval-stochastic  mixed-integer  programming 
(FIMP)  method,  improving  the  robustness  of  the  IDSS  for  real- 
world  applications.  Thus,  it  can  be  used  as  an  efficient  tool  for 
analyzing  and  visualizing  impacts  of  CET  policies,  regional  sus¬ 
tainable  development  strategies,  and  C02  emission  reduction 
measures  in  an  interactive,  flexible  and  dynamic  context. 


6.  Conclusions 

In  this  study,  a  full-infinite  interval-stochastic  mixed-integer 
programming  (FIMP)  method  has  been  developed  for  planning 
carbon  emission  trading  (CET)  under  dual  uncertainties.  FIMP 
incorporates  two-stage  stochastic  programming  (TSP)  and  inter¬ 
val  full-infinite  programming  (FIP)  within  a  general  framework. 
FIMP  has  advantages  in  uncertainty  reflection  and  policy  analysis, 
especially  when  the  input  parameters  are  provided  as  crisp 
intervals,  functional  intervals  and  probabilistic  distributions.  This 
study  attempts  to  introduce  CET  into  the  Beijing’s  electric  power 
system  (EPS),  and  tries  to  give  analysis  of  following  questions: 
(i)  how  to  determine  initial  distribution  of  permits  and  sector 
coverage  for  Beijing’s  domestic  power  market,  (ii)  how  to  quantify 
the  final  system  cost  of  each  power  conversion  technology  under 
CET,  and  (iii)  how  to  ascertain  the  amount  of  electricity  supply 
and  C02  mitigation. 

Then,  a  FIMP-based  carbon  emission  trading  (FIMP-CET)  model 
has  been  formulated  for  planning  C02  emissions  with  trading 
scheme  under  dual  uncertainties  of  electric  power  system  (EPS)  in 
Beijing.  This  is  the  first  attempt  to  introduce  CET  into  Beijing’s  EPS 
to  mitigate  C02  emissions  with  optimization  model.  With  the  aid 
of  model,  the  schemes  about  energy  supply,  electricity  supply, 
carbon  quota,  and  capacity  expansion  have  been  obtained.  Results 
are  helpful  for  supporting  (a)  evaluation  or  improvement  of 
allocation  patterns  of  energy  resources,  (b)  formulation  of  local 
policies  regarding  energy  consumption,  electricity  supply,  eco¬ 
nomic  development  and  energy  structure,  and  (c)  resolving  of 
conflicts  and  interactions  among  C02  reduction,  economic  cost, 
system  reliability  and  energy-supply  security. 

The  results  obtained  suggest  that  the  study  city  should  making 
efforts  to  improve  its  policy  in  the  following  several  manners. 
Firstly,  encourage  renewable  energy  (such  as  hydropower, 
wind  power,  photovoltaic  power  and  biomass  power)  actively, 
such  that  full  implementation  of  these  measures  can  alleviate 
the  contradiction  among  energy  supply,  energy  demand,  and  C02- 
reduction.  Secondly,  improve  C02-reduction  efficiency  (e.g., 
installation  of  C02-reduction  facilities)  and  increase  energy- 
conversion  efficiency,  which  can  evolve  the  international  influ¬ 
ence  of  the  city  fundamentally.  Finally,  the  study  of  C02  capture 
and  storage  technology  should  be  strongly  advocates  for  achiev¬ 
ing  to  the  city’s  commitment  reduction. 
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